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Abstract: We show how to construct lattice sigma models in one, two and four dimen- 
sions which exhibit an exact fermionic symmetry. These models are discretized and twisted 
versions of conventional supersymmetric sigma models with N = 2 supersymmetry. The 
fermionic symmetry corresponds to a scalar BRST charge built from the original super- 
charges. The lattice theories possess local actions and in many cases admit a Wilson term 
to suppress doubles. In the two and four dimensional theories we show that these lattice 
theories are invariant under additional discrete symmetries. We argue that the presence 
of these exact symmetries ensures that no fine tuning is required to achieve N = 2 super- 
symmetry in the continuum limit. As a concrete example we show preliminary numerical 
results from a simulation of the 0(3) supersymmetric sigma model in two dimensions. 
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1. Introduction 

Supersymmetric field theories exhibit many remarkable properties. Chief among these are 
the cancellations which occur between boson and fermion contributions in the perturbative 
calculation of physical quantities. These cancellations eliminate many of the divergences 
typical of quantum field theory and are at the heart of its use to solve the gauge hierar- 
chy problem Q . Additionally, supersymmetric versions of Yang Mills theory, while more 
tractable analytically than their non-supersymmetric counterparts, exhibit many of the 
same features such as confinement and chiral symmetry breaking Q. Finally, supersym- 
metric SU {N) gauge theory in the large N limit has been proposed as a candidate for 
M-theory g. 

Given then, the theoretical and phenomenological interest in supersymmetric theories, 
it is perhaps not surprising that a good deal of effort has gone into attempts to study 
such theories on spacetime lattices see, for example, [Q, |5| and the recent reviews by Feo 
and Kaplan j^, 0. However, until recently these efforts mostly met with only limited 
success. The reasons for this are well known - generic discretizations of supersymmetric 
field theories break supersymmetry at the classical level leading to the appearance of a 
plethora of relevant susy breaking counterterms in the effective action. The couplings to 
all these terms must then be fine tuned as the lattice spacing is reduced in order that the 
theory approach a supersymmetric continuum limit. This problem is particularly acute in 
theories with extended supersymmetry which contain scalar fields. A notable exception to 
this scenario is = 1 super Yang Mills in four dimensions in which it was shown that only 
one coupling need be tuned to restore supersymmetry in the continuum limit H . Because 



of this a good deal of numerical work has been devoted to this model ^ . In |11, IC] we 
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showed that in certain cases this fine tuning problem could be evaded by constructing a 
lattice theory which retained an element of supersymmetry. A similar philosophy has been 
followed in recent papers by Kaplan et al. |12, 13, 14|. The latter work allows constructions 
of lattice super Yang Mills models by orbifolding a large N matrix model. In |15] we showed 
that the models treated in namely supersymmetric quantum mechanics and the 

complex Wess-Zumino model, were actually examples of a wider class of supersymmetric 
theory - namely theories which could be twisted to expose a scalar fermionic symmetry. 
This fermionic symmetry can be written in terms of a charge Q, formed from combinations 
of the original supercharges, and obeying the BRST-like algebra = 0. We will refer to 
it as a twisted supersymmetry. Furthermore, the resulting twisted actions generically take 
the form S = 5A where 5 is the variation generated by Q and A is some function of the 
fields. They resemble pure gauge fixing terms. It is important to realize that in fiat space 
this twisting operation can be considered as merely a relabeling of the fields of the theory 
and hence does not alter its physical content. 

This twisting procedure is very well known in the literature, as the first stage in the 
construction of topological field theories (see for example the reviews and |17] and 



references therein). Indeed Witten pioneered the technique when he twisted N = 2 super 
Yang Mills to arrive at the very first topological field theory |18|. The final transition 
from twisted supersymmetric theory to true topological field theory necessitates restricting 
the set of physical observables to only those invariant under the BRST symmetry. In the 
language of the parent supersymmetric field theory this constitutes a projection to the 
ground state of the theory and is responsible for ensuring that the topological field theory 
does not contain local degrees of freedom. In our constructions of lattice supersymmetric 
models we will not impose this restriction and our actions are simply discrete versions of 
the continuum twisted supersymmetric actions. 

As we have remarked, the structure of these twisted actions suggests that the associated 
topological field theories can be obtained by a gauge fixing procedure. This is indeed 
the case and we will construct our lattice actions by writing down an appropriate BRST 
symmetry and making a suitable choice of gauge fixing condition. We will show that the 
resulting lattice models enjoy a number of nice properties. They are local, yield a = 2 
supersymmetric model in the classical continuum limit and most importantly maintain an 
exact (twisted) supersymmetry. We argue that the presence of this exact symmetry should 
protect the theory from the usual radiative corrections which plague lattice supersymmetric 
models. The construction is sufficiently general that we can, with certain restrictions on 
the target space, formulate the models in one, two or four dimensions. These restrictions 
on the target space are nothing more than the usual restrictions required for the theories 
to exhibit extended supersymmetry. 

In the next section we describe the procedure used to construct these lattice sigma 
models focusing initially on the case of a one dimensional base space (yielding a gener- 
alization of the quantum mechanics model considered earlier fl^). The following section 
shows how to generalize the model to two dimensions leading to a lattice version of the 
usual supersymmetric sigma models with extended supersymmetry. In this case the re- 
quirement that the BRST invariant model possess an interpretation as a twisted form of 
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a supersymmetric theory requires that the target space be Kahler. We then discuss in 
some detail perhaps the simplest example of such a model ~ the 0(3) sigma model. Some 
preliminary numerical results which support our reasoning are then shown. Finally, the 
natural generalization to four dimensions is given which requires that the target manifold 
be hyperKahler. 

Most of the results derived in this paper are not new but have existed in the literature 
of topological field theory for a decade or more. The crucial observation we make here 
is that many of them can be taken over without significant modification to the lattice. 
The purpose of this paper is to explain how this comes about in as pedagogical a way as 
possible. Further discussion of continuum topological field theories and their relationship 
to supersymmetric theories can be found in the excellent review [^] and references therein. 

2. A BRST invariant model on a curved target space 

Consider a set of real bosonic fields (p^{cr) corresponding to coordinates on an A^-dimensional 
target space equipped with Riemannian metric gij{(f)). The coordinates a parametrize a 
base space whose dimension and other properties will be determined later. In addition 
to let us introduce real anticommuting fields iI^^{(t) and r]i{a) and a further bosonic 
(commuting) field Bi{a). These fields will transform as vectors in the target space. Let 
us further postulate the following set of transformations parametrized by a single (global) 
grassmann parameter ^ 

5<j)' = ^V* 
6il>' = 

5m = dB,-r,,Tli;^ 



The quantities F*^ and i?*^; are the usual Riemannian connection and curvature. It is clear 
that = (5^^/^* = 0. It is straightforward, (the details are given in appendix A) to prove 
the additional results 

= 

5^Bi = (2.2) 

(2.3) 



Thus the symmetry is nilpotent. Indeed, as was shown in |1£], these transformations can 
be regarded as a BRST symmetry which arises from quantizing a purely bosonic model 
with trivial classical action 5*01(0) = 0. This classical action is invariant under local shifts 
in the field (j)^ — > + e*. BRST quantization is then employed to define the quantum 
theory. Thus, in this approach, the anticommuting fields if)"^ and r/j are to be interpreted 
as a ghost and antighost field respectively while the field Bi is a Lagrange multiplier field 
used to implement a (yet to be specified) gauge fixing condition. Notice that none of these 
properties depend on the nature of the base space. 
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Having chosen a field content and an associated BRST symmetry it now remains to 
choose an action. Again, the gauge fixing perspective suggests the fohowing natural choice 

= a6 7?, (^iV' (<^) _ ^g'^B,^ ) (2.4) 

The vector N^{<f>) corresponds to an arbitrary gauge fixing condition on the bosonic field 
Notice that the exact nature of the integral over the base space need not be specified 
at this point. The coupling constant a corresponds in this language to a gauge parameter. 
Expectation values of observables which are invariant under the (5-symmetry will not depend 
on the value of a. Notice also that only the nilpotent property of the 6 operator is needed 
to show the the action is BRST invariant. Carrying out the variation leads to the following 
expression for an invariant action 

S = aj^ I^BiN' - ^g'^BiBj - r^iVkN'^j^ + -^Rjimk^f rf^^^^^ (2.5) 

where the symbol V indicates a target space covariant derivative. Integrating out the 
multiplier field B yields the on shell action 



{]^gi,N'N^ - ri.VkN'^'' + (2.6) 



This action is manifestly general coordinate invariant with respect to the target space 
coordinates. It is also invariant under the on shell fermionic transformations 

5^' = Cip' 
5il:' = 

H = e [m,N^ - rijTl^i^'') (2.7) 

Up to this point we have left the choice of the target space vector field arbitrary. In the case 
of a one dimensional base space with coordinate a there is a natural choice for A^* ((/>) = 
This choice ensures that the bosonic action will be quadratic in derivatives which is a 
minimum requirement if this model is ultimately to be interpreted as a supersymmetric 
theory. In this case the action reads 

S = aj da Qg,^.^^ - r^,2_^^ + ii?,,„fcr/'r?V"^') (2.8) 



where 



is the pullback of the target space covariant derivative. This action corresponds to super- 
symmetric quantum mechanics (in Euclidean time) in which the fields take their values on 
some non-trivial target space. It is a straightforward generalization of the model consid- 



ered earlier [10|. It is clear that the special choice of gauge condition we made in making 
contact with the supersymmetric theory corresponds to finding the Nicolai map for the 
supersymmetric theory [pO|]. 
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Perhaps the most important feature of this continuum construction of the (twisted) 
supersymmetric model is that it can be transcribed trivially to the lattice without mutilat- 
ing the fermionic symmetry - simply by replacing the continuum derivative in iV' = ^ 
by a suitable finite difference operator defined on a one dimensional lattice in which the 
continuum coordinate a is replaced by the discrete index t where a = ta with t = 1 ... AT 
and a is the lattice spacing. 

^ ^A+4l, = i {6t>,t+a - St',t) 4 (2-10) 

In this expression we employ a forward lattice difference operator to ensure no doubles 
appear in the spectrum of bosonic states. Notice that the (twisted) supersymmetry then 
guarantees that no such states can appear in the fermion spectrum either. This statement 
is certainly true classically, however we believe that it is also remains true when quantum 
fluctuations are taken into account. Replacing continuum integrals with sums over lattice 
points X in the usual manner we are led to a lattice action of the form. 

= (^5i.A+0^A+</^ - mD^' + \Rjimkrfrf^^'il^^^ (2.11) 

where 

D'^' = A+V^ + r^fc {A+(t^) (2.12) 

While this lattice action is invariant under the twisted supersymmetry it is no longer 
generally coordinate invariant in the target space. This derives from the fact that the 
lattice expression for AT* {(p) no longer yields a target space vector when the base space 
derivative is replaced by a finite difference operator. However, the deviation of N'^ from 
a true vector field decreases smoothly to zero as a ^ cx) and the continuum limit is 
approached. Thus in the vicinity of such a continuum fixed point it appears that the 
lattice theory approaches a generally covariant continuum theory possibly perturbed by 
irrelevant operators. We are currently studying this issue in more detail. 

Notice that the identification of this model with supersymmetric quantum mechanics is 
very straightforward in the case of a one dimensional base space - the twisting procedure 
is trivial in this case. It merely requires one to identify the ghost and antighost with 
the physical "fermion" fields ^ , ^ respectively. In higher dimensions the identification is 
somewhat more complicated since the ghost and antighost will turn out to be related to 
different chiral components of the physical fermions. Notice also that to make contact with 
the supersymmetric theory we do not impose the usual physical state conditions associated 
with BRST quantization - namely that the BRST charge annihilate the physical states. 
Thus the model is not topological. Nevertheless, the BRST symmetry, is exact and will 
guarantee the vanishing of a set of associated Ward identities corresponding to the twisted 
supersymmetry. We now turn to the generalization of this model to describe a field theory 
in two dimensions. 
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3. Twisted two dimensional sigma models 



The model we discussed in the last section is sufficiently general that it can be used to 
generate theories defined on a two dimensional base space. In essence all we have to do is 
to replace the "gauge condition" N^{(p) by some suitable function. To lead to a quadratic 
bosonic action it should involve a single derivative of the field (j) with respect to the base 
space coordinates. A simple guess for the continuum form would be 

Here, the vector A^* must pick up an extra index a = 1,2 corresponding to the two possible 
directions in the base space. Notice that the presence of this extra index implies that 
both antighost and multiplier field also pick up an extra index Tji — > f]ia 

and Bi — > Bi^. 

It is easy to verify that the fermionic transformations are unaffected by the addition of 
this extra index as is the nilpotent property of the BRST operator corresponding to those 
transformations. Actually this choice will not do. It is clear that if we are to arrive at a 
twisted supersymmetric model the number of degrees of freedom carried by the antighost 
must match that of the ghost field (in the end they will turn out to correspond to different 
chiral components of the physical fermions). Thus we must require the antighost r/ and 
multiplier field B to satisfy some condition which halves their number of degrees of freedom. 
The natural way to do this is to introduce projection operators p(~) and P^"*") and require 
that T7*° and P*" satisfy certain self-duality conditions 

pi-)r^ = 

pWr/ = 7? (3.2) 



One choice for these projectors is 



prf^ = \{m^ (3-3) 

Here, Jj must be a globally defined tensor field on the target space and is the usual 
antisymmetric matrix with constant coefficients (at this point we shall think of the base 
space as flat). In order that P(+)2 = p(+), p(-)2 = p(-) and p(+)p(-) = p(-)p(+) = the 
tensor must be antisymmetric and square to minus the identity. Manifolds possessing 
such a structure are called almost complex and have even dimension. At this point we 



must be careful to make sure that the BRST transformations we introduced earlier eqn. 2.1 



are compatible with these self-duality conditions given in eqn. 3.2. Explicitly, the BRST 
transformations will need no modification if 

5 [P;fv^^) = H'" (3.4) 

This latter condition is satisfied provided the almost complex structure is covariantly con- 
stant Vyfc Jj = 0. Manifolds possessing a covariantly constant almost complex structure are 
termed Kahler. Using these projectors the continuum gauge condition now becomes 

^ia ^ pi'^ + Q^^j (3,5) 
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Just as in the previous section we choose a a pure gauge-fixing term as action 

S = a5 j (faVh (^7]ia (^AT^" - ^B'"^ j (3.6) 

where the change in the coefficient of the multipher term is simply to generate a conven- 
tional normalization for the boson kinetic term. In this expression we have now allowed 
for a general metric h^/^ in the base space (in which case we must interpret as an almost 
complex structure on the base space). Carrying out the variation and integrating out the 
multiplier field leads to the on shell action 

- Via^kN'^i^'' + ^i^i.-femr/^VaV'V") (3-7) 

In this expression gij denotes, as before, the metric on the target space. Notice that this 
invariance of the action under the BRST symmetry is true independent of the choice of 
either base or target space metric provided they are both Kahler. This observation lies at 
the heart of the use of these models to construct topological field theories. Furthermore, 
for such manifolds, the projector, buried in the ghost-antighost term, may be taken out of 
the covariant derivative to act harmlessly on the self-dual antighost on the left. This term 
then takes the form 

riio^D^^' (3.8) 

where the base space covariant derivative acting on the ghost field arises from the pullback 
of the target space covariant derivative acting on the gauge function. In detail 

= d^^i + T),^d'^(jy>'il}'' (3.9) 

Furthermore, in the case of a Kahler manifold the piece of the boson action involving Jj is a 
topological invariant. The physical interpretation of this model is clearer if we go to complex 
or Kahler coordinates in the target space. The original target space fields (f)^,i = 1 . . . 2N 
are replaced with complex fields c/)^ ,1 = 1 . . . N and their complex conjugates {(f)^)* = (j)^ ■ 
In these coordinates the only non-zero components of the metric take the form Qj-j and 
cfjj = {(jfjY ■ Furthermore, the tensor Jj locally takes the form i5j for the unbarred 
coordinates and for the barred coordinates. The target space interval can be written 

ds'^ = 2gjjd(t>^d(ty^ (3.10) 

Similarly we can adopt complex coordinates to parametrize the base space which we can 
take as flat. The original real coordinates ai, ct2 are replaced with the complex coordinates 

o"± = o"i lb ia2 with h-\ = h |- = 2 and base space line element dl"^ = 2h^ da'^da~ In 

these coordinates the gauge condition becomes simply the Cauchy-Riemann condition 

d+cp^ = (3.11) 
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solutions of which are holomorphic functions. In such complex coordinates the self-duality 
condition on the antighost field rj implies 



riL =0 

7]l = (3.12) 

Similarly the action takes the form 

S = a j d^a (2h+-gjjd+(f)^d-(p^ 

- h+-g,jr]iD_iF -h+-gjjilD+i:^ + h+-R,jjjr,irlt^iF^ (3.13) 

Notice that the only non-zero terms of the Riemann tensor for Kahler manifolds are pre- 
cisely of the form shown and guarantee that this term is real. The covariant derivatives 
acting on the ghost fields are given by 

D+i^J = d+t' + r;^i9+(/<^V'' (3.14) 

Again, for Kahler manifolds parametrized in complex coordinates the only non-zero con- 
nection terms are of the form and its complex conjugate Tjj^- In terms of these 
complex fields the BRST transformations become 

S(j)^ = ^ip^ (3.15) 
(JV'^ = (3.16) 
5r^i = e {2d+^' + T'jKV+i^') (3.17) 

This action and associated fermionic symmetry is in agreement with Witten's original 
construction of topological sigma models restricted to the case of Kahler target spaces [^l| . 

At this point the physical interpretation of the model is not yet transparent. We need 
to rewrite the theory in terms of physical fermion fields - that is we need to untwist the 
theory. The first stage in this process consists of recognizing that the the operator 5+ can 
be thought of as the Weyl operator in Euclidean space. To see this consider the chiral 
representation of the Dirac matrices in Euclidean space 



71 



Thus the free Dirac operator in this basis is just 





Since Weyl spinors are just complex numbers in two dimensions we see that solutions of 
the gauge condition are in one-to-one correspondence with solutions of the Weyl equation 
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for right handed spinors. Indeed we see that the part of the action |3.13| quadratic in the 
ghost-antighost fields can be rewritten in this chiral basis as 



where 



^' -i-f-D^-^Qj-j (3.18) 

With this relabeUng we can see that the model may be identified with the usual supersym- 
metric sigma model with N = 2 supersymmetry with the gauge parameter a reinterpreted 



as the physical coupling a = p- (see |16],|27]). The crucial ingredient which made this 



final identification possible was the requirement that the target manifold be Kahler. This 
allowed the construction of a covariantly constant projection operator which was compat- 
ible with the BRST symmetry and which could be used to define a self-dual gauge fixing 



condition eqn. 3.5. The structure of the latter then gave rise to an equation for the ghost 
fields that could be interpreted as a Dirac equation for a chiral fermion. Turning this ar- 
gument around it is clear that only supersymmetric models with Kahler target spaces can 
be twisted to reveal actions invariant under a BRST symmetry. 

It is clear that the twisting procedure for these two dimensional sigma models can be 
regarded as essentially a decomposition of the original Dirac field into chiral components. 
To clarify this issue consider the N = 2 supersymmetry algebra 

{Qa+,Ql3-] = la/S^fi 

where +/— denote quantum numbers +1/2 and —1/2 of the internal R-symmetry with 
generator R. The original rotation group has a single generator J with quantum numbers 
in a spinor representation a = + and a = —. The supercharges are then denoted Q±±. 
When the theory is twisted we replace the original generator of rotations J hj J + R. It is 

now clear that two supercharges have spin zero Ql = Q-\ and Qr = Q |-. Furthermore, 

it follows from the original superalgebra that these new charges obey the BRST-like algebra 

qI = Qr = {Ql,Qb] = ^ 

We can then use say Ql to construct our twisted supersymmetric action. The twisting 
procedure does not always yield such a trivial change of variables. Consider, for example, 
four dimensional super Yang Mills theories in which the fields transform as representations 
of SU{2)l X SU{2)fi X SU{2)j where / labels the isospin or R-symmetry. The twisting 
procedure amounts to constructing a new rotation group SU{2)l x SU{2)j and decompos- 
ing all fields into representations of this group. After this decomposition a scalar nilpotent 
charge is produced corresponding to the trace of where i and a are the original isospin 
and spinor indices. In this case the twisted fermion fields correspond to fields which trans- 
form as a scalar, vector and selfdual tensor under the new rotation group. We refer the 



reader to |18| for more details. 

Returning now to the sigma models, it is trivial to see that the invariance of the action 
under the twisted supersymmetry is retained when I replace the continuum action with an 
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appropriate lattice action. However, unlike the case of quantum mechanics, we cannot do 
this merely by replacing all continuum derivatives by forward difference operators - it is 
easy to see that the kernel of the (free) 2D lattice Dirac operator constructed this way still 
contains extra states which have no continuum interpretation^. Instead we must proceed 
in a way similar to the complex Wess-Zumino model [11| and introduce an explicit Wilson 



mass term. Of course it is not obvious that the addition of such a term is compatible with 
the topological symmetry in the case of a curved target space. However, in ||2^, ^ it was 
shown that indeed it is possible to add potential terms to the twisted sigma models while 
maintaining the topological symmetry. In complex coordinates the possible terms are 

AS = X^V^Vi + AV^V/FjV'^ - hi+~r]\ViV-jrii (3.20) 

Here, is a holomorphic Killing vector and an arbitrary parameter. A Wilson term 
would correspond to the choice = imw(t>^ where 

miy = ^(A+Ar + A+A,-) (3.21) 

Many Kahler manifolds possess such a holomorphic Killing vector (for example the CP^ 
models considered in the next section). Thus for these models the doubled modes can be 
removed without spoiling the Q-exactness of the lattice action. 

In detail, the transcription for latticization now entails replacing the continuum deriva- 
tive in the gauge condition by a symmetric finite difference operator 

nI_ = A^cj)^ (3.22) 

where A^ = Af + iAf and A'^ = i (A+ + A^) together with the appropriate potential 
terms following from eqn. 3.2C with = imw(t>^ ■ 



The potential terms should eliminate the doubles from both fermion and boson sectors 
while preserving lattice rotational invariance of the theory at any lattice spacing. Thus far 
we have left the constant free. In the topological field theory language it corresponds to 
a gauge parameter. However, if we untwist the model to a theory of Dirac fermions there 
is a natural choice which preserves the Lorentz invariance of the model A'^ = h'^ = ^ for 
a flat base space. Using this and the explicit form of the Killing vector allows us to rewrite 
the additional piece in the action as 

2 

-AS" = imw4>^) imw4>i) + '^limw'^^ (3.23) 

a 

Where the decomposition of the Dirac fleld in terms of the ghost and antighost was given 
in eqn. 



3.19 



Finally, we turn to a discussion of the renormalization properties of these two dimen- 
sional models. Speciflcally we would like to know whether this lattice model approaches the 
continuum theory with the full N = 2 supersymmetry automatically as the lattice spacing 
a — > 0. The requirement that the action is invariant under the twisted supersymmetry 



^We thank Joel Giedt and Erich Poppitz for this observation 



-10- 



ensures that any local counterterm induced via quantum effects must take the form of a 
gauge fixing term^. Furthermore, simple power counting arguments lead us to conclude 
that the only (marginally) relevant terms take the form 

O = S {v^'^fi, ((/.) {dc.^ + B^)) (3.24) 

where we have reverted to the original formulation involving real fields. Notice that we 
have written this operator in continuum language assuming a restoration of Poincare in- 
variance in the base space. We have also assumed that the latter symmetry additionally 
enforces general covariance in the target space as we discussed in the one dimensional case. 
General covariance ensures that fij is a tensor which may then be taken to represent a 
quantum renormalization of the target space metric tensor. This counterterm structure 
would be consistent with a lattice model which exhibits = 1 supersymmetry in the con- 
tinuum limit. The restoration of full N = 2 supersymmetry appears to require additional 
constraints. Luckily, such constraints are present in the form of additional discrete symme- 



tries of the lattice action. Consider the classical action in Kahler form given in eqn. 3.13 
It is trivial to see that this action is also invariant under the transformations 

TjL -irjL (3.25) 

Actually, this additional symmetry arises from the Kahler structure of the target space 
appearing in the classical action. In pH] it was shown that the N = 2 supersymmetric 



action was invariant under the transformations •^j^"' where is the Dirac spinor 

introduced earlier in connection to the untwisting of the model. The tensor Jj is just the 
(covariantly constant) complex structure characteristic of a Kahler manifold. After twisting 
and adopting complex coordinates we obtain the symmetry shown in eqn. |3^ . This 



additional symmetry of the lattice model then ensures that only counterterms compatible 
with a Kahler target space survive in the quantum effective action. But as was shown 



in [25 1 any model with = 1 supersymmetry and a Kahler target space automatically 
possesses N = 2 supersymmetry. Thus we expect that no additional fine tuning is needed 
to regain the full supersymmetry of the continuum model. 

In the above argument we have only considered radiatively induced operators which 
take the form 50. One might worry that other invariant operators which are not of this 
form could be produced via quantum effects. In fact the list of such operators for the 2d 



a model is very short - it consists of just one element (see |21]) 

O = iJ^j^'i;^ 



^Actually, for the model with potential terms corresponds to a target space coordinate transformation 
along the Killing vector. This is not an exact invariance of the lattice action but yields a supersymmetry 
breaking term which vanishes like a'^. In a superrenormalizable theory this term cannot generate relevant 
supersymmetry breaking terms and the above analysis is still essentially correct 
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In the continuum limit (where the Wilson term vanishes) the chiral symmetry of the bare 
Lagrangian prohibits such a term and hence we can safely neglect it in our analysis. 

We turn now to an explicit example of these ideas - the 0(3) supersymmetric sigma 
model. 



4. Lattice 0(3) nonlinear sigma model 



The CP^ models are perhaps the simplest example of sigma models possessing Kahler 



target spaces. They have been extensively studied in the literature ||25|, in part 

because of their similarities to gauge models in four dimensions. These models contain 
N complex scalar fields acting as coordinates on a complex manifold whose metric is 
locally derived from the Kahler potential F = In ( 1 + u^u^ ) in the usual way 



1 9' --'jj 



The case = 1 is especially interesting as it corresponds to the usual 0(3) supersymmetric 



sigma model [27]. In this case the metric, connection and curvature are easily verified to 
be 

9uu = ^ (4.2) 
217 

= = -- (4.3) 

RuUUU 9uudu^uU A (4'4) 

where p = 1 + uu. In this case the supersymmetric lattice action including Wilson terms 
takes the form 



S = 



X c 1 1 c — 1 c* 

-^A^uA^u + —{mwu){mwu) - -—^r]D_ip - -—^f]D_^_ip 



(4.5) 



where a factor of two has been absorbed into the coupling a and we have simplified our 
notation by replacing r/_|_ — > and r/_ rj. The explicit form of the lattice covariant 
derivative is _ 

Dl = Al-j {Alu) (4.6) 

To proceed further it is convenient to introduce an auxiliary field a to remove the quartic 
fermion term (this field a has nothing to do with the original base space coordinates which 
have been replaced by integer lattice coordinates x). Explicitly we employ the identity 



J I)ae~"(^"""^^'"^+^'''^) = e^'"''^'^ (4.7) 
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where is the number of lattice sites. Thus the partition function of the lattice model 
can be cast in the form 

Z = J DuDaDriDije-'^^''^''^'^''^^ (4.8) 
where the action is now given by 



S = 



(4.9) 



"1 1 1 

p"^ p'^ 2 

1 - 1 1 - 1 

and the covariant derivative is modified to include a coupling to the auxiliary field a 

D+ = Al- — {Alu) + a (4.10) 

Notice that the factor of coming from the introduction of the auxiliary field is canceled 
by a similar factor resulting from the original gaussian integration over the multiplier 
field B. Thus, this partition function should be independent of the coupling constant a 
which plays the role of a gauge parameter from the perspective of the BRST construction. 
We will present numerical results confirming the a independence of Z shortly. Using the 




decomposition eqn. p.l9| the Dirac operator M now takes the form 

-(ziJ+ji imw I 

To simulate this model we have to reproduce the fermion determinent arising after inte- 
grating out the anticommuting fields. This leads to an effective action of the form 

S = aSsiu, a) - ^Tr In (^a'^M^u, a)M{u, a)) (4.12) 

where 5s (it, a) denotes the local bosonic pieces of the action and we have shown the de- 
pendence on coupling a explicitly. This representation of the fermion determinent requires 
that the latter be positive. This is certainly true in the continuum limit when M is positive 
definite and we observe it to be true in all our simulations at finite lattice spacing. For small 
lattices we have employed a HMC algorithm to yield an exact simulation of this non-local 
effective action |2^ . This algorithm requires a full inversion of the fermion matrix at every 
time step and is prohibitively expensive for large lattices where we have instead employed 
both a stochastic Langevin scheme and the so-called R-algorithm to simulate the system 
pO| ], |pl|]. In these cases physical quantities exhibit systematic errors 0(At)^. Typically, 
for the data presented here we have used A = 0.05 which yields a systematic error smaller 
than our statistical error. 

A stringent test of the a-independence of Z is gotten by measuring the expectation 
value of Sb- It should be clear that 

.'2^^,^^Ss)-?^ (4.13) 
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where V denotes the number of lattice sites. The results are given in the following table 
which shows data for ^ < Sb > from runs at a variety of couplings a for lattices of 
size 4x6, and 8x8. The runs comprise 1.7 x 10^ trajectories and 1 x 10^ trajectories 
respectively. 

Notice that while these numbers are statistically differ- 
ent from one at small a they rapidly approach unity as a 
increases. This is in agreement with our expectation that 
the Q-symmetry is only broken by terms which vanish as a 
power of the lattice spacing. Most importantly no fine tun- 
ing is needed to see a restoration of this symmetry in the 
continuum limit a — ^ oo. There is another way to under- 
stand this coupling constant independence of the partition 
function. In a conventional supersymmetric theory the par- Table 1: "^2V^'' 

tition function gives a representation of the Witten index 
W. The latter can be written as 

W = Tr{-lfe-'^^ 

where F is the fermion number operator. The topological character of W then relies on 
the bose-fermi pairing of all states with non negative energy. Thus the a independence 
of Z, guaranteed by the Q-cxactness of the twisted action, appears to imply that the 
Hamiltonian of the twisted system will indeed possess the exact degeneracy required of a 
supersymmetric theory. 

As a further check on the supersymmetry of the model we have studied the simplest 
two point functions involving boson and fermion fields. The boson correlator G^{t) = 
Re (< u{t)u{Q) >) projected to zero spatial momentum is shown in figure 1. as a function 
of the (Euclidean) time coordinate t for an 8 x 8 lattice at a = 2.0. The curve shows the 
result of a simple fit to the functional form a-|-6cosh (m^ {t — T'/2)) yielding = 0.417(7) 
as an estimate of the lowest lying boson mass state (T denotes the length of the time axis) . 
In this and the following fits we exclude the t = data as it will contain contributions 
from higher mass states. The fermion correlator projected to zero spatial momentum 
G[j{t) =< "^j(t)^'j(0) > should take the form 

Gfj{t) = A{t - T/2)Sij + iB{t - T/2)eij 

where A and B are even and odd functions and we have taken the Dirac gamma matrix 
in the time direction to correspond to a2- All of our data very accurately reproduces this 
spinor structure. Figure 2. shows a plot of ImG^]^(t) as a function of t again for an 8 x 8 
lattice at a = 2.0. The solid curve shows a simple fit to the form B = a sinh (m^^ (t — T/2)) 
with rriQ-^ = 0.41(1). Similarly the function ReG^o(*) the same run is shown in figure 
3. together with a simple fit to a hyperbolic cosine yielding a consistent estimate for the 
fermion mass m^Q = 0.409(1). Thus the fits are consistent with the equality = as 
required by supersymmetry. 

We arc currently extending these simulations to larger a and bigger lattices to study 
the approach to the continuum limit. Notice that a great deal is known theoretically about 



a 


4x6 


8x8 


0.5 


0.953(5) 




1.0 


0.969(4) 




1.5 


0.987(4) 




2.0 


0.999(4) 


0.985(3) 


2.5 


1.002(4) 




3.0 


1.007(5) 


0.995(3) 
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Figure 1: G^{t) for a = 2, L = 8 
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Figure 2: G^^{t) for a = 2, L = 8 



the mass spectrum of this model in the continuum |32] and it would be very interesting 
to check the results of our simulations against this work. Furthermore, it is possible to 
write down lattice Ward identities for both the exact and broken supersymmetries. Study 
of these should give information on the question of the restoration of the full N = 2 
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supersymmetry in the limit of vanishing lattice spacing. These questions are also under 
investigation. 

5. Lattice Sigma Models in Four Dimensions 

Let us turn now to four dimensions. Again, the appropriate gauge condition (cp) = 
must take the form of a (base space) derivative of the target space bosonic field (p'^. To 
ensure that the corresponding antighost r/jQ, and multiplier field Bia possess the correct 
numbers of degrees of freedom this gauge condition should employ some projector to enforce 
a self-duality condition. In four dimensions this condition should reduce the number of 
degrees of freedom by a factor of four. In two dimensions the projector depended on an 
almost complex structure for both base and target spaces both of which were required to 
be Kahler. The local complex structure associated with a Kahler manifold then ensured 
that the ghost-antighost action could be reinterpreted as action for Dirac fermions in a 
chiral representation. In the same way a local almost quaternionic structure is necessary to 
allow for a representation of the four dimensional Dirac gamma matrix algebra. Manifolds 
with an almost quaternionic structure possess three independent tensors (J^Yj o, = 1 ... 3 
satisfying the quaternionic algebra 



jajb ^ _^abj ^ ^abcjc 



(5.1) 



In terms of these tensors the projector we need looks like [22] 




(5.2) 
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Here, the tensors J" are almost complex structures in the target space and j"" are the 
corresponding base space quantities. The gauge condition in the continuum now reads 

^ia ^ pia+Qp^ (5.3) 

and as before the antighost and multiplier fields are required to be self-dual under the 
action of this projector 

P+r? = ri (5.4) 

In a way similar to two dimensions the requirement that the self-duality conditions are 
compatible with the BRST symmetry requires all three almost complex structures to be 
covariantly constant. Four dimensional manifolds admitting three such covariantly constant 
tensors are called hyperKahler. 

It is straightforward to find an action corresponding to fixing this gauge condition 

^S = asJ d^aVhmcc [n"^ - ^i?'") (5.5) 

Carrying out the variation and integrating out the multiplier field (taking into account the 
self-duality condition) yields the following on shell action 

S = aj d^aVh (^h^^gijda<l>%(t>> + \ {fT^ (J"),,- 5«0^a^<^' 

- h-f'gijriiDf.^P^ + ^%fei/i"^j,r;^^VV') (5-6) 

Notice again the presence of the base space covariant derivative in the ghost-antighost 
term. For hyperKahler manifolds the term in the bosonic action containing the tensors 
J" and j'^ is again a topological invariant. It is possible to untwist this model to produce 
a Lorentz invariant theory of Dirac fermions in a way completely analogous to the two 
dimensional case. To make this explicit it is convenient to introduce a notation that 
naturally emphasizes the quaternionic nature of the hyperKahler space. 
First, group the original AN target fields in sets of four 



9 ,9 ,9 ,9 ,9 ,9 ,9 ,9 ,-,9 ,9 ,9 ,9 (5.7) 

^',1=1 ^1,1=2 ^i,I=N 

then the gauge condition condition, Pj^^c?^^ = can be rewritten (locally) as 

\ d^^i - d,^l - ei^^afidc,^'^ = ^ ■ ^ 

($^ really means ^jt^ 1 = 1, ...,N and fx = 1,2,3,4). To see that this is true 

let us choose the following representation for the J"'s (likewise for the j"'s) 

Ji=7=f^^^ ° ] J2 = j=f ° J3 = i^=('0 (5.9) 

\ ia2 Wi / Us i ^ ' 
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Equations (|5.q ) are equivalent to the Cauchy-Fueter equations that generahze to four di- 



mensions the Cauchy-Riemann equations of two dimensions |22]. Finally let us adopt 
quaternionic or hyper-Kahler coordinates by defining 



= - - + $1 and d+ = Idi + Jd2 + Kd^ + 84 

then equations (|5.8| ) for the gauge condition take the simple form, 

d:<^^ = 



(5.10) 



(5.11) 



Solutions to this equation define so-called triholomorphic functions. The action for the 4D 
model may then be written in a way analogous to the two dimensional case eqn. |3.13| by 
allowing all fields to take values in the quaternions and recognizing that a single quater- 
nionic antighost (and its conjugate) t]^ = + Irjf^ + Jr]2^ + Krj^^ survives the self-dual 
projection. 

The physical interpretation of this model necessitates untwisting. The procedure mim- 
ics the 2D case. First, it is easy to show that the solutions to the 4D gauge condition 
correspond to solutions of the Weyl equation for right handed spinors just as in the two 



dimensional case |22|. Indeed, if we take the chiral Euclidean representation of the four 
dimensional Dirac matrices, 



(5.12) 




yiai j 

the free Dirac operator again becomes simply 

7.9 = 





(5.13) 



(where d- = —Idi — J 82 — Kd^ + d^) and we adopt a common representation of the 
quaternionic structure in terms of the Pauli matrices / = —idi, J = —icr2 and K = —ia^. 
Then the action of this operator on a right handed Weyl spinor. 



/ \ 





(5.14) 



with, 



= -i- 



(5.15) 



again just yields the set of equations (|5.8|) . Thus, the local quaternionic structure of the 
hyperKahler target space automatically allows us to encode the Clifford algebra of the 
Dirac matrices and ensures that the ghost and antighost field correspond to the chiral 
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components of a single Dirac fermion. Thus the untwisted model corresponds to the usual 
N = 2 theory based on hyper-multiplets [P 

As with the one and two dimensional models it is trivial to discretize this theory 
while preserving the twisted supersymmetry by simply replacing the derivative operator 
appearing the gauge condition eqn. 5.3 by a symmetric difference operator. Such a choice 
will give rise to both bosonic and fermionic doublers. To remove such states from the 
spectrum we again need to add a Wilson term in the form of a simple mass operator to the 
theory. By analogy with the two-dimensional case it would seem that this would only be 
possible if the target space of the model admits a triholomorphic Killing vector linear in 
the field. Several examples of such theories are discussed in [28| for hyperKahler manifolds. 
For such manifolds we can add a Wilson type operator of the form m]y(f) while maintaining 
the Q-exactness of the lattice action. 

Following the arguments given earlier the exact supersymmetry will help protect the 
theory against radiative corrections. Indeed, as in two dimensions, the 4D lattice model 
will possess three additional discrete symmetries corresponding to acting on the spinors 
of the theory by each of the three independent complex structures. In the continuum 
these additional discrete symmetries generate additional supersymmetries. Clearly the 
requirement that the lattice effective action respect these additional symmetries will force 
the target manifold to remain hyperKahler under renormalization. This in turn should 
ensure the full N = 2 supersymmetry is achieved at large correlation length without further 
fine tuning. 
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6. Conclusions 



In this paper we have shown how sigma models in one, two and four dimensions with N = 2 
supersymmetry may be discretized in a way which preserves an exact super symmetry. The 
argument rehes on the well known twisting procedure used to produce topological field 
theories in the continuum. Here, this argument is turned around and topological construc- 
tions are used to generate twisted supersymmetric theories which may then be discretized 
while preserving an exact (twisted) supersymmetry. The twisted supersymmetry takes the 
form of a BRST symmetry and the resulting twisted supersymmetric action resembles a 
gauge fixing term. Notice that while the method we use to construct our lattice actions 
parallels continuum topological field theory constructions our theories are not topological. 
They are to be regarded as simply relabeled forms of conventional supersymmetric actions. 
Indeed, we have untwisted the models explicitly to reveal the physical fermion interpre- 
tation. The ghost and anti-ghost fields just correspond to different chiral components of 
a Dirac fermion field. Notice that from the point of view of simulation this untwisting 
procedure is somewhat unnecessary - the twisted fields generate precisely the same fermion 
determinant as the physical fields and it is only this object that we must replicate in any 
simulation. 

The models we construct are local and, at least for target spaces possessing suitable 
holomorphic Killing vectors, can be rendered free of doubling problems while preserving 
the Q-exactness of the lattice action. Our preliminary numerical results for the CP^ 
model lend support to these claims. Furthermore, it seems plausible (though we have no 
proof) that the restoration of Poincare invariance leads to a restoration of target space 
reparametrization invariance in these theories. If this is so then we have argued that the 
exact supersymmetry together with certain discrete symmetries should automatically lead 
to a restoration of full N = 2 supersymmetry in the continuum limit without additional 
fine tuning. We are currently investigating these issues in more detail. 

Of course one of the most important questions is whether these ideas may be extended 
to theories with a gauge symmetry. In the continuum it is indeed possible to construct 
twisted versions of a variety of super Yang Mills model (the original example being afforded 
by = 2 theory in four dimensions [^]). However, there is an important distinction 
between these twisted gauge models and both the sigma models described here and the 
Wess-Zumino models that were considered earlier |15|. The relevant BRST symmetry is 
now nilpotent only up to gauge transformations - so if we want to latticize such models it 
appears we must be careful to preserve gauge invariance. While this manuscript was being 
prepared we received a preprint in which significant progress in this direction had been 
achieved. 
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A. Appendix 

In this appendix we want to show that 



indeed from (Ejl 



S\ = 

6^Bi = (A.l) 



but, 



where {'0'",^'^} = has been used. Thus, 

= -l^xvjR'uki^'i^' + l^xvjR'iiki^'^' 

= (A.3) 

now consider 

= e \x{BiT',mi'"' - IvuR'^jimi^'i^nnui^' + B.d^vlx^^i^' 



2 



1 



= ^XB,{TI^T\, + d^T\^i,^^^ - -ixB,R\,^^'^'' 

^\iXvAdmR\,^ + R\,^V\ + Vi^R\^)^^^'^^ (A.4) 
using (A.2), the two first terms cancel each other. The last term can be written as, 

(where we used the anti-symmetry of the Grassman's), 
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(as the two added terms vanish when contracted with 

= (A.6) 
(by Bianchi identity). Thus 5^ = i.e. is a nilpotent operation on the multiplet (^, 77, B, ip). 
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